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ABSTRACT 

We searched for quasars that are ~ 3 mag fainter than the SDSS quasars in 
the redshift range 3.7 < z < 4.7 in the COSMOS field to constrain the faint end 
of the quasar luminosity function. Using optical photometric data, we selected 
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31 quasar candidates with 22 < i' < 24 at z ~ 4. We obtained optical spectra for 
most of these candidates using FOCAS on the Subaru telescope, and identified 
8 low-luminosity quasars at z ~ 4. In order to derive the quasar luminosity 
function (QLF) based on our spectroscopic follow-up campaign, we estimated 
the photometric completeness of our quasar survey through detailed Monte Carlo 
simulations. Our QLF at z ~ 4 has a much shallower faint-end slope (/? = 
— 1.67lo 17) than that obtained by other recent surveys in the same redshift. Our 
result is consistent with the scenario of downsizing evolution of active galactic 
nuclei inferred by recent optical and X-ray quasar surveys at lower redshifts. 

Subject headings: cosmology: observations — quasars: general — surveys 

1. Introduction 

The quasar luminosity function (QLF) is one of the most important tools to constrain 
the evolution of supermassive black holes (SMBHs). A number of luminous quasars have 
been found up to z ~ 6 by the Sloan Digital Sky Survey (SDSS; Fan et al. 2000, 2001, 2003, 
2004, 2006; Richards et al. 2006; Goto 2006; Jiang et al. 2008, 2009) and 2dF quasar surveys 
(Boyle et al. 2000; Croom et al. 2004). However, low-luminosity quasar^] at z > 4 have not 
yet been well studied. In contrast, the QLF at z < 3 is well quantified in a wide luminosity 
range (e.g., Croom et al. 2009) and is best represented by a double power law (e.g., Boyle 
et al. 1988; Pei 1995). Accordingly the bright-end slope a has been already well measured; 
the typical measured value is around —3.1 at z < 2.4, and flattens to a > —2.37 by z = 5 
(Richards et al. 2006). More interestingly, recent studies on the optical QLF show that the 
activity in low-luminosity active galactic nuclei (AGNs) peaks at a lower redshift than that 
of more luminous AGNs (Croom et al. 2009). By assuming that the brighter AGNs have the 
more massive SMBHs, this process can be interpreted as AGN (or SMBH) downsizing. The 
AGN downsizing is also reported by X-ray surveys (Ueda et al. 2003; Hasinger et al. 2005; 
see also Brusa et al. 2009). However, the lack of low- luminosity quasars at z > 4 leaves the 
faint-end slope of the z > 4 QLF very poorly constrained. Consequently it is not understood 
how low-luminosity quasars evolve at high redshifts, or if downsizing is also present in the 
earlier universe. 

Motivated by these issues, we have searched for low luminosity quasars at z ~ 4 in the 



1 Although quasars have been defined as relatively high-luminosity AGNs with Mb < — 21.5 + 5 log 
ho historically (Schmidt & Green 1983), we do not adopt such a solid criterion to distinguish these two 
populations. 
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COSMOS field (Scoville et al. 2007). Throughout this paper we use a A cosmology with f2 
= 0.3, fi A = 0.7, and the Hubble constant of H = 70 km s _1 Mpc -1 . 



2. The Sample 
2.1. The Data 

We select the quasar candidates by using color and morphology criteria. We use the 
official COSMOS photometric redshift catalog (Ilbert et al. 2009). This catalog covers an 
area of ~ 2 deg 2 and contains several photometric magnitudes including the total magnitude 
(MAG_AUT0) of the z'-band and 3" diameter aperture magnitudes of the CFHT w*-band as 
well as Subaru Suprime-Cam g'-, r'-, and i'-bands. Details of the Suprime-Cam observations 
and the COSMOS photometric catalog are given in Taniguchi et al. (2007) and in Capak 
et al. (2007), respectively. The 5<x limiting AB magnitudes are u* = 26.5, g' = 26.5, r' 
= 26.6, and i! = 26.1. Since we also use the Advanced Camera for Surveys (ACS) catalog 
(Koekemoer et al. 2007; Leauthaud et al. 2007) to separate galaxies from point sources, our 
survey area is restricted to the area mapped with ACS (1.64 deg 2 ). 

2.2. Quasar Candidate Selection 

A quasar at z ~ 4 shows a Lyman break in its spectral energy distribution (SED) that 
falls between the wavelengths of the g' and r' filters, making a g' — r' color redder than their 
r' — i' color. We utilize this characteristic to select candidates of low-luminosity quasars. 
Here the typical quasar colors as a function of redshift are necessary to define the quasar 
color-selection criteria. Therefore, we created model quasar spectra following the procedure 
generally adopted (e.g., Fan 1999; Richards et al. 2006), and derived the u* — g' , g' — r', and 
r' — %' colors of the model quasars at redshifts between < z < 6. For the model quasar 
spectra we assume typical values of the power-law slope a v (/„ oc v~ av ) and Lya rest-frame 
equivalent width (EW) to be 0.46 and 90A, respectively (Vanden Berk et al. 2001). In Fig. 
1, the simulated colors of the model quasars are shown in the r' — i' versus g' — r' diagram. 
We select our candidates of quasars using the following criteria: (1) 22 < i'(MAG_AUT0) < 24, 
(2) r' - i' < 0.42(o' - r') - 0.22, (3) u* - g' > 2.0, and (4) g' - r' > 1.0. where we use 
MAG_AUT0 measured by SExtractor (Bertin & Arnouts 1996) in eq. (1) instead of 3" aperture 
magnitude, because the target luminosity should be calculated from total magnitudes, not 
from aperture magnitudes. The criterion (2) is used to select quasars efficiently without 
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significant contamination of stars. To remove possible foreground contaminations further, 
we introduce the additional criteria (3) and (4). These latter two color thresholds (i.e., 
criteria (3) and (4)) are determined by taking empirical color distributions of quasars at 
z ~ 4 into account (Richards et al. 2006). We also exclude spatially-extended objects 
defined by Leauthaud et al. (2007). As a result, we obtained 31 quasar candidates among 
7318 point sources with 22 < i'(MAG_AUT0) < 24. 

3. Spectroscopic Observation 

The spectroscopic follow-up observations of the quasar candidates have been carried 
out at the Subaru telescope with FOCAS (Kashikawa et al. 2002) on 7-11 January 2010. 
We used the 300 grating with the S058 filter, whose wavelength coverage is 5800A < A obs 
< 10000A (1040A < A rcst < 1590A at z = 4). We used a O'.'8-width slit, resulting in a 
wavelength resolution of R ~ 700 (Av ~ 430 km s _1 ) as measured by night sky emission 
lines. The typical seeing was ~ Qf!7. We observed 28 of the 31 candidates, and obtained 
useful spectra for 23 of them. The remaining 5 objects show no signal, possibly due to 
insufficient exposure time (the typical exposure time of these objects is 1800 seconds). 

Standard data reduction procedures were performed using IRAF. After sky subtraction, 
we extracted one- dimensional spectra by adopting an aperture size of 1'.'8. The spectra of 23 
objects were flux-calibrated using spectrophotometric standard stars. We found that eight 
show strong and broad Lyct and C IV emission lines, suggesting that these eight objects 
are indeed quasars at z ~ 4. Of the remaining 15 objects, one spectrum shows narrow 
C IV, He II, O ill], and C ill] emission lines (type-II quasar) at z ~ 3.5. The spectra of 
three other objects show only narrow emission lines without any high-ionization lines such 
as C IV present, being consistent with Lya emitting galaxies at z ~ 4. The remaining 11 
objects appear to be consistent with Galactic late-type stars (from K4 to M0 stars). The 
photometric properties and the measured redshifts of the identified quasars are summarized 
in Table 1. An example of the reduced spectra of the identified quasars is shown in Fig. 2. 

4. Quasar Luminosity Function 

Our spectroscopic run found eight quasars at z ~ 4. Accounting for the possibility that 
some of the photometric candidates without spectra could be also quasars (based on the 
success rate of our spectroscopic run), the "corrected" number of the quasars in our survey 
are (4.35 ± 2.01) x 10" 7 Mpc~ 3 mag" 1 for -24.09 < M 1450 < -23.09 and (7.32 ± 2.97) x 
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10~ 7 Mpc~ 3 mag" 1 for -23.09 < M li50 < -22.09, respectively. To derive the QLF based 
on these results, we need to compute the effective comoving volume. The effective comoving 
volume of the survey is calculated as: 

pz=co 

Vrff(wv) = dQ C(mi>,z)—dz, (5) 

Jz=o dz 

where dfl is the solid angle of the survey and C(m,i',z) is the photometric completeness. 
Here the photometric completeness is calculated by modeling quasar spectra as described in 
Section 2.2, but taking also the intrinsic variation of the continuum slope and the equivalent 
width (EW) of emission lines into account. We assume a Gaussian distribution of the power- 
law slope a v (/„ oc v~ av ) and Lya EWs, with means of 0.46 and 90A (the same as those 
in Section 2.2), and a standard deviation of 0.30 and 20A, respectively (Vanden Berk et al. 
2001). We created 1000 quasar spectra at each Az = 0.01 for the redshift range < z < 6. 
The effects of intergalactic absorption by neutral hydrogen were corrected by adopting the 
extinction model of Madau (1995). Then, we calculated the colors of the model quasars 
in the observed frame. We put the simulated quasar images into Subaru FITS images and 
measure their colors. Applying the color-selection criteria (1) - (4), we infer the photometric 
completeness. The redshift range where the completeness is moderately good is 3.7 < z < 
4.7. Specifically, the inferred completeness is ~ 0.7 for quasars with i' = 22.5, and ~ 0.5 for 
those with i' = 23.5. Given the effective comoving volume, the comoving space density and 
its standard deviation are calculated through the standard procedure (see, e.g., Fan et al. 
2001). The obtained QLF is plotted in Fig. 3. Although the redshift range of SDSS data 
(4.0 < z < 4.5; Richards et al. 2006) is slightly different from our study (3.7 < z < 4.7), we 
apply the weighted least-squares fit to the space density of quasars at z ~ 4 inferred by our 
study and SDSS, adopting the following double power-law function: 

(b( M 7 ) = W*45o) (R , 

^1^1450, &) 10 OA( a +l)(M 1450 -M* i50 ) + 10 0.4(/3+l)(M 14 50-^ 450 ) ' ^ 

where a, /3, $*(Mi45o), and Mj* 450 are the bright-end slope, the faint-end slope, the normal- 
ization of the luminosity function, and the characteristic absolute magnitude, respectively. 
Among the four parameters, the bright-end slope (a) is fixed to be a = —2.58 based on 
the z ~ 4.2 SDSS results (Richards et al. 2006). The best-fit parameters are $(MZ 450 ) = 
(3.20 ±0.24) x 10~ 7 Mpc~ 3 mag" 1 , M* 450 = -24.40 ±0.06, and (3 = -1.67lg;#. The fitting 
result is shown in Fig. 3. 

Here we compare our QLF at z ~ 4 with another QLF for the same redshift recently 
reported by Glikman et al. (2010). They searched for quasars at z ~ 4 using the data of 
the NOAO Deep Wide-Field Survey (NDWFS; Jannuzi k Dey 1999) and the Deep Lens 
Survey (DLS; Wittman et al. 2002), whose total area is 3.76 deg 2 . By combining their 
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own data with the SDSS results, they derived a faint-end slope of /3 — —2.3 ± 0.2. This is 
significantly steeper than our result (/5 = — 1.67^° '{],). This discrepancy could be possibly 
caused by the fact that only six spectra have been obtained for a total of 117 photometric 
quasar candidates at R > 23 in the sample of Glikman et al. (2010) and there is a possibility 
that a large fraction of their photometric quasar candidates at R > 23 are contaminants 
(Lya emitting galaxies and Lyman break galaxies), since such high-z star-forming galaxies 
could be also selected through the Lyman-break selection criteria. While we removed such 
star-forming galaxies by using their ACS images, Glikman et al. (2010) used ground-based 
images to distinguish point sources and extended sources. Note that the number density 
of quasars at R > 23 in Glikman et al. (2010) is roughly consistent with that of Lyman 
break galaxies at z ~ 4 with similar luminosity (e.g., Ouchi et al. 2004). The fact that 
the two luminosity functions are roughly consistent suggest a large contaminant fraction in 
their quasar candidates. It should be also mentioned that Cristiani et al. (2004) reported a 
quasar number count at z ~ 3 — 5, which requires a flattened faint-end slope in the QLF. 

We now discuss the evolution of the quasar space density. To compare the quasar 
space density with previous studies, we need to estimate the quasar space density in the 
same magnitude bins. Therefore we re-calculated the space density of quasars with —23.5 < 
^1450 < —22.5. Fig. 4 shows the redshift evolution of the quasar space density for different 
absolute magnitude bins. Although there are a number of low-luminosity quasar surveys at 
z ~ 3 (Wolf et al. 2003; Hunt et al. 2004; Fontanot et al. 2007; Bongiorno et al. 2007), we 
plot only the results of the 2dF-SDSS LRG and Quasar Survey (2SLAQ; Croom et al. 2009), 
the Spitzer Wide-area Infrared Extragalactic Legacy Survey (SWIRE; Siana et al. 2008), and 
SDSS (Richards et al. 2006), in order to avoid data with large statistical errors. Although 
most studies at z < 3 have suggested consistent results (AGN downsizing), the situation 
becomes rather controversial at z > 3. Specifically, the results of this study and by Glikman 
et al. (2010) show completely different pictures. If the result of Glikman et al. (2010) 
is correct, the high number density of low-luminosity quasars (that may have less- massive 
SMBHs) at z ~ 4 would correspond to the "seeds" of luminous quasars (with high-mass 
SMBHs) at z ~ 2 — 3. However, our results do not show any evidence of such a break down 
of the downsizing scenario at z ~ 4, suggesting that numerous seeds of high mass SMBHs 
should exist at even higher redshift (z > 5). Note that the type-II (i.e., obscured AGN) 
fraction is reported to be higher in lower-luminosity samples at higher redshifts (La Franca 
et al. 2005). Therefore one may suspect whether the possible abundant obscured population 
affects the analysis of the (un-obscured) quasar number density evolution, since we are now 
focusing on high redshift low-luminosity AGNs. However, given the redshift and luminosity 
of our sample, the inferred obscuration fraction is ~ 0.3 at most (Hasinger 2008). This 
suggests that our analysis is not significantly affected by the obscured population. To obtain 
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a more conclusive understanding of the quasar evolution in the early universe, extensive 
surveys of high- z low-luminosity quasars are definitely required. The next-generation wide- 
field prime-focus camera for the Subaru Telescope (Hyper Suprime-Cam; Miyazaki et al. 
2006) and Large Synoptic Survey Telescope (LSST; Ivezic et al. 2008) will provide us with 
such opportunities in very near future. 

We would like to thank the Subaru staff for their invaluable help and all members of 
the COSMOS team. TM , KM , and YI are financially supported by the Japan Society for 
the Promotion of Science (JSPS) through the JSPS Research Fellowship. 
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Table 1: Redshift and optical photometry of the eight spectroscopically confirmed z ~ 4 



quasars 





Redshift 


u* a 


g' 




%' 


i! (MAG_AUT0) 


Exp. Time b 


Number 


z c 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(sec) 


1 


3.89 


>27.49 


25.05 


23.56 


23.51 


23.35 


6000 


2 


4.14 


>27.49 


24.75 


23.10 


22.79 


22.55 


6000 


3 


3.56 


>27.49 


25.82 


24.42 


24.09 


23.91 


3600 


4 


4.20 


>27.49 


27.27 


24.09 


23.68 


23.45 


1800 


5 


3.86 


27.42 


24.62 


23.45 


23.21 


23.04 


2400 


6 


3.65 


27.00 


24.14 


22.54 


22.35 


22.16 


3600 


7 


4.45 


>27.49 


25.23 


22.79 


22.22 


22.03 


3600 


8 


4.16 


>27.49 


24.94 


23.27 


22.95 


22.78 


7200 


a When the 


■u*-band 


magnitude 


is fainter than 2a limitin 


g magnitude (= 


27.49), the 2a 



lower limit is given. 

b Total on-source exposure time in the FOCAS spectroscopic observation. 
c z SDec is based on the C IV emission line. 



This preprint was prepared with the A AS IATgX macros v5.2. 
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Fig. 1. — Two color diagram between r' — il and g' — r' used for our quasar selection. 
Green points are point sources with 22 < i'(MAG_AUT0) < 24. Blue triangles (8 objects) are 
spectroscopically confirmed z ~ 4 quasars. The red line is the median track of modeled 
quasar colors as a function of redshift. The black solid line shows our photometric criteria 
used to select quasar candidates. 
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Fig. 2. — Reduced spectrum of the quasar No. 8 (upper panel) and typical sky spectrum 
(lower panel). Dotted lines show the expected wavelengths of quasar emission lines: Lya 
A1216, N v A1240, O I A1304, Si iv A1400, C IV A1549, and He II A1640. 
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Fig. 3. — The z ~ 4 quasar luminosity function. The red squares show our results. Triangles 
are the data of NDWFS+DLS using M 1450 from K-corrected magnitude (blue) and directly 
from the quasar spectra (magenta) (Glikman et al. 2010). Green squares are the SDSS 
results for the z = 4.25 bin (Richards et al. 2006). The red line shows the fit to our results 
(i.e., COSMOS) and the SDSS points. The black dash-dotted line shows the fit to the data 
from the NDWFS+DLS plus SDSS study (Glikman et al. 2010). 
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Fig. 4. — The redshift evolution of the quasar space density. Red, green, blue, magenta, 
light blue, orange, black, purple lines are M 1450 = —20, —21, —22, —23, —24, —25, —26, and 
-27, respectively. Dotted lines show the combined 2SLAQ, SDSS, SWIRE, NDWFS and 
DLS QLF. Dashed line shows the combined this work (i.e., COSMOS) and 2SLAQ QLF. 



